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ABSTRACT 

We present results on the thermal and kinetic Sunyaev-Zel'dovich (SZ) effects from a sequence of high 
resolution hydrodynamic simulations of structure formation, including cooling, feedback and metal in- 
jection. These simulations represent a self-consistent thermal model which incorporates ideas from the 
'prc-hcating' scenario while preserving good agreement with the low density IGM at z ~ 3 probed by the 
Ly-a forest. Four simulations were performed, at two different resolutions with and without radiative 
effects and star formation. The long-wavelength modes in each simulation were the same, so that we can 
compare the results on an object by object basis. We demonstrate that our simulations are converged 
to the sub-arcminute level. The effect of the additional physics is to suppress the mean Comptonization 
parameter by 20% and to suppress the angular power spectrum of fluctuations by just under a factor 
of two in this model while leaving the source counts and properties relatively unchanged. We quantify 
how non-Gaussianity in the SZ maps increases the sample variance over the standard result for Gaussian 
fluctuations. We identify a large scatter in the Y-M relation which will be important in searches for 
clusters using the SZ effect (s). 

Subject headings: cosmic microwave background - cosmology: theory - galaxies: clusters: general - 
large-scale structure of universe - methods: numerical 



1. INTRODUCTION 

On angular scales below 10', the microwave sky carries 
the imprint of large-scale structure in the low-z universe. 
In particular, cosmic microwave background (CMB) pho- 
tons propagating through the universe are inverse Comp- 
ton or Doppler scattered by hot electrons along their path, 
either in dense structures such as clusters of galaxies or 
more generally in hot gas in the intergalactic medium. In- 
verse Compton scattering conserves the number of pho- 
tons but preferentially increases their energy, leading to 
a spectral distortion whose amplitude is proportional to 
the product of the electron temperature and density (or 
pressure). Doppler scattering induces an intensity fluc- 
tuation with the same spectral shape as the CMB itself. 
These effect s were first described by Sunyaev & Zel'dovich 
(|1972| , |1980| ) and are known as the thermal and kinetic SZ 
effects, respecti vely ( for recent reviews see Rephaeli 1995 
and Birkinshaw 1999). 

The thermal effect is one of the primary sources of sec- 
ondary anisotropies in the CMB on small angular scales. 
The change in the (thermodynamic) temperature of the 
CMB resulting from scattering off non-relativistic elec- 
trons is 
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where x — his/kTcMB — ^/56.85 GHz is the dimension- 
less frequency, and the second expression is valid in the 
Rayleigh- Jeans limit which we shall assume henceforth. 
The quantity y is known as the Comptonization parame- 



ter and is given by 
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where the integral is performed along the photon path. 
Since T e Tqmb the integrand is proportional to the 
integrated electron pressure along the line of sight. 

The kinetic SZ effect arises from the motion of ionized 
gas with respect to the rest- frame of the CMB. The result- 
ing temperature fluctuation is AT/T = — 6, where 
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measures the magnitude of the effect along the line of 
sight if v r is the radial peculiar velocity of the gas (pos- 
itive if the cluster is moving away from us, negative if it 
is moving towards us). The different dependence on fre- 
quency of the two effects can, in principle, be used to dis- 
entangle them observationally, though the primary CMB 
anisotropies provide "noise" for the kinetic SZ signal. 

In this paper, we report results of a sequence of high res- 
olution hydrodynamic simulations designed to study the 
SZ signal on small angular scales. The imprint of the SZ 
effects on the CMB has be en st udied by a numb er of au- 
thors (Ostr iker & Vishniac 1986| ; Persi ct al. 1995; da Silva 
et al. |200C| , |2001| Refre gier et al. |2000| ; Seljak et al. pOOO ; 
Kay, LiddlclT Thorn as |200l| ; Refregier & Teissier (2001; 
da Silva et al. |200l|; Zhang, Pen & Wang |2002|). Ours 



are the first such simulations to incorporate the effects of 
radiation and star formation, including gas cooling and 
feedback from supernovae and galactic winds. As such 
they extend our earlier work on th e subject, reported in 
Springel, White & Hernquist (|200l|). 
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The outline of this paper is as follows. In Section 0, we 
briefly describe our simulations and the techniques we use 
to compute SZ maps. In Section ||, we present our results 
for the power spectra of thermal and kinetic SZ effects and 
discuss the effects of non-Gaussianity on the sample vari- 
ance. We show SZ source counts as a function of source 
strength and source size and how SZ source properties cor- 
relate with the intrinsic properties of the cluster. Finally, 
we summarize and discuss our results in Section ^. 

2. METHOD 

We chose as our cosmological model the currently fa- 
vored ACDM cosmology. In particular, we adopted f2 m = 
0.3, n A = 0.7, H = lOO/ikms^Mpc- 1 with h = 0.7, 
l^B = 0.04, n — 1 and a% = 0.9. This model yields a rea- 
sonable fit to the current suite of cosmological constraints 
and as such provides a good framework for making realis- 
tic predictions. The simulation domain was a 100 ft, _1 Mpc 
cube with periodic boundary conditions. We ran smoothed 
particle hydrodynamic (SPH) simulations at two different 
mass resolutions corresponding to 2 x 144 3 and 2 x 216 3 
particles, with initial conditions constructed such that the 
large-scale modes were equal for both resolutions, allow- 
ing us to perform a resolution study on an object by ob- 
ject basis without being affected by cosmic variance. For 
each mass resolution, we performed a simulation that com- 
puted only adiabatic]] gas physics and shock heating, and 
one that also followed the radiative cooling and heating 
processes of the gas in the presence of a UV background 
radiation field together with star formation and feedback 
processes in the dense collapsed gas. In Table [y we list 
some basic numerical parameters of the four simulations 
that we carried out. 

Radiative cooling and heat ing was modeled as in Katz, 
Weinberg & Hernquist ( 1996| ), with a photo-ionizi ng flu x 
similar to the one advocated by Haardt & Madau ( 1996| ). 
Radiation from quasars reionizes the Universe at redshift 
z ~ 6 in this scena rio (fo r details, see e.g. Dave, Hernquist, 
Katz & Weinberg 199S| ) . We have treated star formation 
in the framework of an effective multi-phase model for the 
interst ellar m edium (ISM) , as detailed by Springel & Hern- 
quist ( [2002a ). In this model, rapidly cooling gas of high 
overdensity is assumed to give rise to the formation of cold 
clouds, embedded in an ambient hot medium. The cloud 
material forms the reservoir of baryons available for star 
formation. Massive stars are taken to explode as super- 
novae on a short timescale, and to release their energy as 
heat to the ambient medium of the ISM. Supernovae are 
also assumed to evaporate cold clouds, thereby establish- 
ing a tight self-regulation cycle for the star-forming ISM. 

In addition, we have included a phenomenological model 
for galactic outflows in our simulations, motivated by the 
large body of observational evidence for t he ubi quitous ex- 
istence of such galactic winds (e.g. Martin 199E). Galactic- 
scale outflows are thought to play a crucial role in the 
enrichment of the low-density intergalactic medium with 
heavy elements and for the global regulation of star for- 
mation, and they may have influenced the properties of 
intragroup and intracluster gas in a decisive way. Note 

1 ln what follows, we will refer to the runs that ignored radiative 
effects and star formation as 'adiabatic' but it should be kept in 
mind that these simulations did include shock heating and so are 
not strictly speaking adiabatic. 



that without the inclusion of such strong feedback pro- 
cesses, simulations tend to substantially overpredi ct th e 
luminosity density of the universe (Balogh et al. 2001 ). 
However, the model for galactic winds that we use here 
is able to reduce the global efficiency of cooling and star 
formation su fficient ly to match these constraints (Springel 
& Hernquist ^002^ ). 

In this model, each star-forming galaxy drives a wind 
with a mass-outflow rate equal to two times its star for- 
mation rate, and with a wind- velocity of v w = 484kms _1 . 
These choices are deliberately 'extreme' in the sense that 
the total kinetic energy of the wind is of the order of the 
energy released by the supernovae. However, these param- 
eters are quite typical for the observed properties of out- 
flows from star-forming disks. Whether or not the wind 
can escape from a galaxy depends primarily on the depth 
of its dark matter potential well. Halos with virial tem- 
peratures below ~ 10 6 K can lose some of their baryons to 
an outflowing wind, while more massive halos will contain 
the wind, such that the wind becomes progressively more 
unimportant for the regulation of star formation in more 
massive objects. 

The simulations considered here are part of a larger sim- 
ulation program which shows that the particular feedback 
model we employ allows a numerically converged predic- 
tion for the total star formation density of the universe as 
a function of epoch, provided all star-forming halos above 
virial temperatures of 10 4 K can be resolved. Of course, 
the simulations analyzed here cannot reach such a high 
mass resolution, because they need to model a volume 
large enough to contain rich clusters of galaxies. It is thus 
inevitable that part of the star formation is lost due to res- 
olution limitations. This can be a sizable effect. Springel 



& Hernquist (2002c) show that about 50% of all stars are 
born in halos with virial masses below 10 12 /i _1 M0, which 
corresponds roughly to the mass of the smallest halos in 
our G4 simulation with reliable estimates for the star for- 
mation rate. Note that the 'unresolved' star formation in 
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Fig. 1. — The evolution of the mean (mass weighted) temperature 
in the 4 simulations. Parameters for the simulations can be found 
in Table hi an 'a' after the name indicates that feedback and cooling 
were not included. 
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small halos is also the one that would predominantly be 
able to drive truly outflowing winds. It should therefore 
be kept in mind that the effects of galactic outflows on the 
intracluster medium are likely going to be underestimated 
by the present simulations. 

The simulations were performed on an Athlon-MP clus- 
ter at the Center for Parallel Astrophysical Computing 
(CPAC) at the Harvard-Smithsonian Center for Astro- 
physics. We used a modified ver sion o f the GADGET- 
code (Springel, Yoshida & White |200l| ), and integrated 
the entropy as the in depend ent thermodynamic variable 
(Springel & Hernquist p002b| ). 

In all our runs, the full simulation box was output ev- 
ery 100/i _1 Mpc along the line-of-sight between redshifts 
z ~ 19 and z — 0. This corresponds to 77 dumps per 
model. We made maps^ in a manner similar to da Silva 
et al. (2000 ) as described in our earlier paper (Springel et 
al. |2001| ). To recap briefly we stacked the various outputs 
from the boxes back along the line of sight, with each box 
randomly translated and oriented in the direction of either 
the x, y or z axis. Then, a grid of 512 2 rays subtending a 
constant angle of 1° from the observer was traced through 
the boxes starting at z = 19. 

We produced maps of the y-parameter, the Doppler 
6-parameter, the projected gas density and the thermal 
brcmsstrahlung emission. In the adiabatic simulations we 
do not track the ionization fraction, x e , dynamically. In- 
stead, we set it to 1.158, corresponding to full ionization, 
for temperatures above 10 4 K and take the gas at lower 
temperature to be neutral. The extra physics models self- 
consistently track x e and this is used in the computation. 

It is also convenient to have a 'cluster catalog' from the 
simulations. We produced this by r unnin g a friends-of- 
friends group finder (e.g. Davis et al. 1985| ) with a linking 
length of b = 0.15 (in units of the mean inter-particle 
spacing). The FoF algorithm partitions the particles into 
equivalence classes by linking together all particle pairs 
separated by less than a distance b. We define the center 
of each halo as the position of the potential minimum, cal- 
culating the potential using only the particles in the FoF 
group. This proved to be more robust than using the cen- 
ter of mass, as the potential minimum coincided closely 
with the density maximum for all but the most disturbed 
clusters. Additionally, the position of the center was very 
insensitive to the presence or absence of the particles near 

2 The missin g fact or of h in our code which caused the maps from 
Springel et al. (2001) to be in error has been corrected here. 



Model N p 
G3 2 x 144 3 
G4 2 x 216 3 



M DM 
2.5 
0.7 



M gas 
0.4 
0.1 



€ 

18 
12 



Table 1 

The parameters of the simulations run. The 
simulation box was 100/l _1 mpc on a side. particle 

MASSES ARE GIVEN IN UNITS OF 10 10 hT 1 Mq. GRAVITY WAS 
SOFTENED WITH A SPLINE ON A COMOVING SCALE e (iN 
/l _1 kpc), WITH THE FORCE BEING FULLY NEWTONIAN 
BEYOND 2.8 e AND THE POTENTIAL OF A SINGLE PARTICLE 
BEING —Gm/e AT ZERO LAG, EQUIVALENT TO A 
PLUMMER-SOFTENING ON THE SAME SCALE. 



the outskirts of the halo, and thus to the precise linking 
length used. For each of the halos we computed the mass, 
-^200 j enclosed within a radius, t^oo, interior to which the 
density contrast was 200 times the critical density. We 
shall sometimes refer to this as the virial radius, and de- 
note its angular extent at the distance of the cluster by 

#200- 

3. RESULTS 

The evolution of the mean (mass weighted) temperature 
in the simulations is shown in Fig. [l]. The formation of 
structure in the universe leads to a steady heating with the 
mean temperature today being approximately 1 /3keV. We 
can see the effects of cooling most clearly at early times 
in the suppression of (T) in runs G3 and G4 with respect 
to G3a and G4a, but by z = most of the differences 
have disappeared. The slight jump in (T) at z ~ 6 is due 
to reionization by quasars in this model, which raises the 
background radiation field and leads to increased photo- 
heating of the gas in the universe. 

We show a typical example of our thermal SZ maps in 
Fig. 0. These maps were produced by coadding a large 
number of partial maps, each giving the contribution of 
one of the boxes that we stacked along the photons' path. 
We made 15 maps for each of the 4 simulations, using the 
same set of random offsets and orientations for each of the 
4 simulations. As we noted before (Springel et al. [2001 ) 
filaments, while prominent in the partial maps, are largely 
hidden in the high level of background arising from the 
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Fig. 2. — One of our Compton y maps, from the simulation G4 
including cooling and feedback as described in the text. The corre- 
sponding map from the adiabatic simulation looks almost identical 
to the eye. 
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sum mation over many of these structures (see also Croft 
et al. |2001| ). 

The mean Comptonization we predict from our 15 maps 
is 2.4 x 10~ 6 and 2.5 x 1CP 6 for the lower and higher 
resolution adiabatic models respectively. This becomes 
2.2 x 1CP 6 and 2.1 x 10~ 6 for the runs with cooling and 
feedback. Thus we see that the effect of cooling is more 
important than heating in aggregate for the gas contribut- 
ing to the mean y distortion and that our m odels easily 
pass the constraint from FIRAS (Fixsen et al. 1996). 



3.1. Power spectra 

For each map we computed a number of statistics, and 
we in general averaged the results for 15 random lines of 
sight to reduce field-to-field variance. Here we present the 
angular power spectrum of the thermal and kinetic SZ 
effects (Fig. |). 

A comparison of the low and high resolution simula- 
tions indicates that for the adiabatic case we have con- 
verged in the angular power spectrum out to i ~ 20, 000. 
These results are also in excellent agreemen t with our ear- 
lier estimates (Springel, White & Hernquist 2001) once we 
take into account the slight differences in the cosmological 
model. For the runs including cooling we are more sensi- 
tive to the effects of finite mass and spatial resolution. The 
higher resolution simulation shows a lower power spectrum 
than the lower resolution spectrum, indicating that cooling 
has been able to proceed more efficiently when the simula- 
tion can correctly resolve denser structures. However, the 
final results are quite insensitive to extra physics included 
in these models, as might have been expected if the bulk 
of the signal is coming from gas outside the core of the 
cluster, and the gas content of the intracluster medium is 
not significantly depleted by gas cooling. In t his re spect 
our results differ from those of da Silva et al. (2001) who 
found a significant effect on the SZ signal if they include 
cooling in their simulations. The addition of star forma- 
tion and energy injection in our simulations stabilizes the 
runaway cooling which otherwise occurs and significantly 
reduces the impact of the extra physics on the SZ effect. 
In the absence of such feedback the cooling is only miti- 
gated by numerical resolution, leading to an excess of cold 
gas and stars in conflict with observations. 

The power spectra we predict are in accord with recent 
measurements of arcminute ang ular scale an isotropy in the 
CMB at BIMA (Dawson et al. [200 1| ; ^002[ ) though lower 
than the de tectio n quoted by CBI (Mason et al. [2002 



Bond et al. 2002] ) from their deep field (the data from 



mosaic observations, Pearson et al. 2002, do not signifi- 
cantly constrain the signal at high-/J! To bring the adi- 
abatic model into agreement with the CBI points would 
require us to e.g. increase <j% to 1-1.2 (or d ramat ically in- 
crease JTgj see e.g. Seljak, Burwell & Pen 200C or Bond 
et al. 2003) . The inclusion of cooling and feedback lowers 
our predictions by ~ 30% on CBI scales, which could be 
compensated for by an additional ~ 4% change in ag. 

We have also estimated the amount by which the evident 
non-Gaussian nature of the SZ signal increases the field- 
to-field variance in the power spectrum. For a Gaussian 
signal the fractional error per mode in the power spectrum 
is v2 ; so we expect the standard deviation in the binned 
power to be y/2/N Ce where N is the number of modes in 
the bin. If we compare this expectation to the standard 
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Fig. 3. — (Top) The angular power spectra for the thermal and 
kinetic SZ effects averaged over 15 maps for each of our 4 sim 
lations. Parameters for the simulations can be found in Table 
an 'a' after the name indicates that feedback and cooling were no" 
included. The upper solid curve shows the primary anisotropy sig- 
nal for this model (including the effects of lensing). The symbols 
the mean SZ signal and the error bars the standard deviation of 
the power from the 15 maps. A dashed line connects the results of 
simulation G4. The lower solid curve is an estimate of the effects 
of patchy reionization assumed to occur at z = 6. The points with 
error bars are observational results from BIMA and CBI (see text). 
(Middle) The fractional error in the thermal effect, per binned C'i , 
obtained from the variance of the 15 simulated fields (open symbols) 
compared to the expectations for a Gaussian sky (filled symbols). 
Symbol types denote two different binnings of the power spectrum. 
The non-Gaussian nature of the SZ signal increases the field-to-field 
scatter over the Gaussian expectation. (Bottom) The angular power 
spectra for 5 of our fields, showing that the ficld-to-field fluctuation 
is primarily in the amplitude, and not the shape, of the spectra. 
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deviation computed from the 15 maps (Fig. [}]), we find 
that the latter is higher by a factor of 2 (at £ ~ 10 3 ) to 
5 (at £ ~ 10 4 ) due to the contributions from the 3- and 



4-point functions (e.g. Meiksin & White 1999, §2). Since 



the non-Gaussianity leads to mode-coupling, the exact ra- 
tio depends slightly on the binning chosen - we show two 
choices. We caution the reader that all of the maps are 
created from the same underlying simulation volume, and 
so are not fully independent. This makes the estimates 
at lower £ somewhat uncertain. We do not have enough 
maps to accurately calculate the off-diagonal elements of 
the correlation matrix (which converge slowly compared 
to the diagonal elements), however the correlations are 
strong. As for the non-linear mass power spectrum, the 
primary fluctuation from field to field is in the amplitude, 
and not the shape, of the power spectrum. We show an 
example from 5 of our fields in Fig. |[ 

Finally it is of interest to ask where the majority of the 
thermal SZ signal is coming from in the angular power 
spectrum - i.e. diffuse gas or gas within halos. To answer 
this we remade two of our maps excluding all gas with den- 
sity smaller than 100 p gas . Gas in halos should have higher 
density than this. The resulting maps are indistinguish- 
able by eye, even 'blinking' between them. There is a small 
(25%) suppression of power for £ < 2000 in the maps with 
the density cut, but the suppression is < 2% by £ = 6000 
where the signal peaks. This suggests (da Silva et al. 2001) 
that even the large-angle signal is coming primarily from 
relatively high density gas along the line-of-sight — the 
virialized halos - and not from the diffuse gas or filamen- 
tary structures (the 'cosmic web'). If we increase the cut 
even further, to 10 3 p gas , we start to significantly reduce 
the thermal SZ signal up to £ ~ 10 4 , indicating that the 
dense cores of the halos do not dominate the thermal SZ 
signal. 

3.2. Source counts 

A glance at Fig. ^ indicates that the SZ effect is dom- 
inated by discrete sources and that the pixel distribution 
is highly non-Gaussian. Thus a power spectral analy- 
sis contains only a small part of the information in the 
map. To obtain an understanding of the SZ source sizes 
and strengths we have used the sour ce detection software 
SExtractor (Bertin & Arnouts 1996) to count individual 
sources, and to measure their strength. We used default 
settings for the detection algorithm of the image process- 
ing software and let it estimate background and 'noise' 
levels automatically. Note that our maps are in principle 
noise-free, so what is interpreted as noise in this procedure 
effectively arises from limitations due to source confusion, 
which leads to non-detections of some of the faintest and 
smallest sources. In a typical analysis with SExtractor, we 
resolve slightly more than half of the total SZ signal into 
discrete sources. 

We measure the strength of a source as the monochro- 
matic brightness change 



S v = I AB V dQ = f(x)B v [ y{9) dfl 
Jn Jn 
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Fig. 4. — Upper panel: the cumulative distribution of source 
fluxes at v = 30GHz. Lower panel: the cumulative distribution of 
source sizes. 



of the CMB integrated over the solid angle of the source. 
Here B v is the Planck spectrum of the primary CMB, and 



6 



f(x) is the spectral function 



i 



e x - 1 



-4 



(6) 



with a; = hv/kTcMB and we have taken v = 30GHz. 

In Fig. ||, we show the cumulative source counts per 
square degree found in the simulations. Again we find re- 
markably good agreement between all of the simulations, 
with the cooling runs predicting fewer sources at the bright 
end by 20-30% and at the faint end by 30-50% (at fixed 
flux). This is not terribly surprising at the bright end 
where the source counts are more massive clusters. In this 
case the bulk of the ICM is at too low a density to cool ef- 
ficiently. However our calculations suggest that this holds 
true even for the relatively weak SZ sources. Cooling and 
feedback thus appear to be a small effect on the clustering 
and number counts of sources, with the gas dominating 
the SZ signal being more prone to cooling than heating in 
our model. 

Each of these bright sources is associated with at least 
one massive halo in the simulation. It is instructive to pur- 
sue this association further. In Fig. ^ we show a scatter- 
plot of source strength vs. halo mass, both intrinsic and 
from the maps. For the 'intrinsic' effect we calculated the 
SZ signal coming from the most massive clusters in our 
box at z = 0. Specifically we summed the SZ signal within 
half of T2oo to obtain the total decrement intrinsic to the 
cluster: 



of redshifts. Second, clusters are not spherical, leading to 
differences in signal when viewed from different orienta- 
tions. Finally, the signal we see is the projection of SZ 
signal from the entire line-of-sight ^ Since the signal is not 
dimmed by distance, objects all along the line-of-sight can 
contribute a significant amount to the signal in any pixel. 
A similar 'confusion' effect operates for opti cal an d weak 
lensing cluster searches (Whit e & Kochanek 2002; White, 
van Waerbeke & Mackey |00|) . Since clusters form in high 
density regions, with a lot of dense material nearby, there 
is a non-negligible chance that some lines of sight will be 
enhanced by non-cluster material. The extreme outliers 
in Fig. are all clusters closer to the observer than about 
600ft.~^4pc, whose virial radii subtend a large angle and 
are most likely to suffer from this projection effect. These 
clusters are marked with crosses in the figure. 

To isolate this effect further we have computed the inte- 
grated SZ signal coming from spheres of gas, centered on 
rich clusters, viewed from many random directions. Specif- 
ically, we first extract all of the gas particles in the simu- 
lation within a distance r cut of the minimum of the clus- 
ter potential. From this set we compute, using Eq. (^), 
the total Y for all particles whose 2D projected radius 
R < T200 /2 for several lines-of-sight. We also compute Y 

4 This contribution was not included in the work of Holder et 
al. (2000) who simulated only gas within the virial regions of their 



y int = ( 7 -i)(i 



*r x \ -> nij u.j 



where 7 = 5/3, Y p = 24%, ot is the Thomson cross sec- 
tion, m e the electron mass, c the speed of light, fij is the 
mean molecular weight of particle j, Xj is the fractional 
ionization (relative to hydrogen), m,j is the mass and u 



the internal energy. The sum is over all particles withi: 
r 20o/2 of the potential minimum in the group, which is the 
region dominating the SZ emission. Since y is dimension- 
less we quote Y mt in units of (/i -1 Mpc) 2 . As the total sig- 
nal scales as M x T and for self-similar clusters T ~ M 2 / 3 
we expect Y"i„t ~ M 5 / 3 which is a good fit to our simula- 
tions. The upper panel of Fig. shows that there is a very 
tight relation between Y lnt and M200, as has been noted 
before by e.g. Holder et al. fl200Cp . It is this relation which 
is behind the claims that a y-selected sample will be close 
to mass limited. 

However, Yj n t is not an observable. To mimic some- 
thing which is closer to what an observation could probe 
we first found all of the clusters in the simulation which 
fell in the line-of-sight of a particular map. For each of 
these we integrated the Compton y parameter values in 
pixels which lay within half #200 of the cluster center to 
obtain Y map = j dfly. As the lower panel of Fig. || shows, 



there is a large amount of scatter in this Y m 



M rela- 



tion, even for our noise free, high resolution maps. This 
scatter is approximately constant as we vary the aperture 
used to define Y map and can be traced to 3 prime sources. 
First, there is evolution in the M — T relation with red- 
shift, and the clusters we have plotted are from a range 

3 The degree of correlation is almost independent of the exact 
aperture used. 
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Fig. 5. — The relation between the total decrement, Y, and cluster 
mass A/200- The upper panel shows the decrement calculated from 
the material within half the virial radius r2oo for each cluster in 
the simulation volume at z = 0. The straight line indicates Yi nt <x 



5/3 
200 • 



The lower panel indicates the decrement obtained directly 
from the noise free maps by summing pixels centered on the known 
locations of clusters (excluding clusters near the edge of the map). 
Open circles indicate clusters whose disk overlapped that of a more 
massive cluster in the map. Filled circles are 'isolated' clusters. 
Crosses indicate clusters within 600h _1 Mpc of the observer. The 
source of the scatter is discussed in the text. 
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from the subset of the particles contained within our FoF 
group. 

Our results are shown in Fig. ^ for r cu t — 5/i _1 Mpc 
and r cut = 10 ft, _1 Mpc. For several clusters there is a large 
scatter between lincs-of-sight depending on the presence or 
absence of material near the cluster. The scatter is larger 
for the larger radius, indicating that it is not just neigh- 
boring material which drives the scatter. If we include 
only those particles within the FoF group we still see a 
large scatter, though the SZ signal is noticeably depressed. 
Thus, some of the scatter is coming from asphericity of the 
cluster, but a non-negligible fraction of the signal is coming 
from non-cluster material along some lines-of-sight. 

The large scatter that we see in the SZ signal at fixed 
mass has implications for SZ searches for galaxy clusters - 
any method which uses primarily SZ flux in target selec- 
tion will not be mass limited. In fact, the steeply falling 
mass functions predicted by hierarchical theories would 
imply that some high-flux candidates are likely to be low 
mass groups with abnormally high SZ signal rather than 
the intrinsically rarer rich clusters. If this is true, a blind 
SZ search should turn up many groups in addition to the 
desired cluster sample. 



including the extra physics on the 1- and 2-point func- 
tions are slight, indicating that the SZ effect is not overly 
sensitive to the (uncertain) thermal history of the gas. 

The maps are noticeably non-Gaussian and dominated 
by discrete sources. This non-Gaussianity manifests itself 
in increased scatter, by a factor of 2-6, in Cg from ficld-to- 
field compared to the Gaussian expectation. 

The properties of the sources in the map are quite in- 
sensitive to the extra physics we have included also, and 
we are thus able to make robust statements about their 
detection. The brightest SZ sources have a density of ~ 1 
per square degree and their detection requires an angular 
resolution of ~ 1', and the capability to separate bright- 
ness fluctuations of size ~ 10 mJy from the primary CMB 
fluctuations. 

We have found that there is significant scatter in the SZ 
fluxes within apertures centered on clusters. This scatter 
has 3 sources: evolution of the mass temperature relation, 
asphericity in the matter distribution and line-of-sight pro- 
jection. The existence of this scatter means that surveys 
using only SZ flux will select objects with a wide range of 
masses, including lower mass groups. 



4. CONCLUSIONS 

We have performed a sequence of high resolution hydro- 
dynamic simulations of structure formation in a ACDM 
model to investigate the effect of extra physics on the ther- 
mal and kinetic Sunyaev-Zel'dovich effects. Including only 
adiabatic gas physics our simulations of the thermal effect 
are converged down to sub-arcminute scales and agree well 
with our earlier work. For the runs with extra physics we 
see an increased efficiency of cooling in the higher reso- 
lution simulation which leads to a lowering of the mean 
Compton y parameter and a slight depression in the SZ 
angular power spectrum. Overall however the effects of 
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Fig. 6. — The scatter in the Yj nt — M relation coming from ma- 
terial close to the clusters. Y\ nt is calculated for each cluster in the 
simulation volume at z = when viewed down 64 random lines 
of sight. Only material within i?2oo/2 in projection was included. 
Crosses include all material with r C ut = 10/i _1 Mpc, open squares 
rcut = 5/i Mpc and dots only the material belonging to the FoF 
group. Points have been offset (horizontally) for clarity. 
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